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Abstract 

The expected sensitivity of Belle is reported for the Dark Photon, A, and Dark Higgs, h! searches, for mass ranges, 
respectively of 0.27 < < 3 GeV/c 2 and 0.54 < my < 10 GeV/c 2 . The Dark Photon and Dark Higgs will be 
searched for in the Higgs-strahlung channels: <? + e~ — > Ah', with h' — > AA and A — > l + l~ (with I — e or //). Results will 
be unblinded soon. 
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1. Introduction 

Ordinary matter represents 4 % of the total energy bud- 
get of the universe. The remainder is believed to be par- 
titioned as either Dark Energy (73%) or Dark Matter 
(23%), but the nature of those components is unknown. 
However, experimental results from direct Dark Matter 
searches, (e.g. DAMA/LIBRA) UHl and other exper- 
imental observations that may be interpreted as devia- 
tions from the Standard Model (e.g. g-2 J5J), can be 
explained in the Dark Matter (DM) models by the inclu- 
sion of an additional interaction, a dark U(l) interaction 
i3- H^l . This interaction, which is mediated by a dark 
U(l) boson, also known as the "Dark Photon", typically 
has very small couplings to Standard Model particles. 

Dark gauge bosons are postulated to have low mass; 
of order MeV to GeV due to astrophysical constraints 
1 17. [HI]. These astrophysical observations include: 
excesses in the cosmic-ray flux of electrons and/or 
positrons above expected background beyond normal 
astrophysical processes and the expected flux of protons 
and/or anti-protons. Dark matter could be charged un- 
der the dark U( 1 ) symmetry group and then the observed 
excess might correspond to dark matter annihilating into 
a Dark Photon A, which in turn decays into l + T (with / 
= e or fj. or possibly r if energetically allowed). 

The ideal tools to discover such particles are there- 
fore not the highest energy hadron collider experiments, 
but lower-energy lepton high-luminosity collider exper- 



iments such as Belle/Bellell and BaBar/SuperB, or ded- 
icated fixed target experiments, several of which are 
planned or already under construction at JLAB (New- 
port News, USA) or at MAMI (Mainz, Germany), for 
example. In Belle, work on dark gauge boson searches 
was started only recently, and has focused on the strate- 
gies proposed by Hill. These proceedings will fo- 
cus on the so-called Higgs-strahlung channel, e + e~ — » 
Ah' and in particular the decay mode with 3e + 3e~ and 
3ji/ + 3/i~ final states. The dark U(l) symmetry group 
could be spontaneously broken, often by a Higgs mech- 
anism, adding a dark Higgs h' (or dark Higgses) to 
these models. The dark photon A can decay into ei- 
ther l + r, hadrons or invisible particles and h' into AA, 
l + r, hadrons or invisible particles. The decay mode of 
the A and h! depends of the mass difference between A 
and h! [21]: (a) to/,- < m&: h' — > invisible particles, (b) 
ma < my < 2mA'. h' — > l + T or hadrons, (c) my > 2m a'- 
h' — > AA. Case (c) will be discussed in this proceedings. 

The A and h' do not necessarily have prompt decays 
I 2ll I24H . The decay length of the dark photon is a 
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function of the dark photon coupling strength to Stan- 
dard Model fermions and is proportional to the inverse 
of the square of the dark photon coupling. The decay 
length of the dark Higgs varies for my > 2m^ between 
being prompt and one meter. In the Higgs-strahlung 
channel, two couplings are involved: the electromag- 
netic coupling of the dark photon to the Standard Model 
fermions, a'; and the dark photon coupling to the dark 
Higgs, au- KLOE and BaBar reported their searches 
on the dark photon and the dark Higgs |25] and ll26ll 
(see also contribution of B. Oberhof to these proceed- 
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ings B27I0 . KLOE focused their search on my < yiia and 
BaBar on my > 2mA (and with A and h' prompt), but 
no signal was found. Only BaBar could set limits for 
the highest mass range, 0.8 < /«/,< < 10.0 GeV/c 2 and 
0.25 < m A < 3.0 GeV/c 2 . BaBar looked at two types of 
final states, fully reconstructed {3l + l~, 2l + 2l~n + n~ and 
2n + 2n~l + l~ where / = e,/i); and partially reconstructed 
(2fi + 2fi~X and e + e~p + [i~X, where X denotes any final 
state other than a pair of pions or leptons). The reason to 
look for partially reconstructed final states is that above 
m A > 1 GeV/c 2 , final states with hadrons dominate as- 
suming BF(A -> e+O + BF(A -> + BF(A -> 
hadrons) = 1 and that ^^^y = R, where R is the 

o~(e + e~ —>hadrons) 



hadronic cross section ratio, R 



BaBar 



cr(e + e — >fJ + fi ) 

reported six candidate events detected: one Apln, two 
47r2/i, two 4n2e and one 4-pX events, in a ~ 500 fb-1 
data set. There were no candidates in the final states 
with six leptons. The number of events detected were 
consistent with background expectations. For low mass 
range of the dark photon: 1 < mA < 300 MeV/c 2 , 
part of the space, coupling strength to Standard Model 
fermions versus dark photon mass, have been excluded 
by Beam dump experiments ( I28I434IO by looking at re- 
action pp — > AX for example. 



2. Experimental setup 

The Belle detector is a large-solid-angle magnetic 
spectrometer, which consists of a silicon vertex detec- 
tor (SVD), a central drift chamber (CDD), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation coun- 
ters (TOF), and an electromagnetic calorimeter (ECL) 
of CsI(Tl) crystals located inside a super-conducting 
solenoid that provides a 1.5 T magnetic field. An iron 
flux-return (KLM) located outside the coil is optimized 
to detect K° L mesons and to identify muons. A de- 
tailed description can be found in 13511. Belle is cur- 
rently being upgraded to Belle II, an upgraded detector 
for operation at SuperKEKB, which will have 40 times 
higher luminosity than KEKB Q. The KEKB collider 
B37I1 . located in Tsukuba, Japan, is the world's highest- 
luminosity electron-position collider. KEKB has pro- 
duced more than one ab~' of data at center-of-mass en- 
ergies corresponding to the T(1S) to Y(5S) resonances, 
and in the nearby continuum. 



3. Particle and reaction identification 

The sub-systems used to identify the electrons are 
primarily the ECL, which measured the energy and 



the CDC, which measured the momentum 13 al . The 
muons are identified by using the KLM and an analy- 
sis combining the measurements of penetration depth, 
the charged track, and the muon cluster matching [39]. 
Events with six charged tracks with three pairs of oppo- 
site charges are considered for analysis. Furthermore, in 
order to maintain a high detection efficiency we require 
that at least three out of the six tracks be identified as 
leptons. 



4. Analysis strategy 

Events with six lepton final states from e + e~ — > 
Ah' — » AAA — » 3/ + 3T (I — e or fj) are reconstructed. 
Energy and momentum conservation is required. The 
invariant mass for each combination of leptons is re- 
quired to be consistent with three distinct A — » l + t . 
Combinations with three "equal" masses (ml, m^ and 
w?,) and mim > 2mu are kept. The "equality" is de- 



fined as follows: m™ ea " 



■3.o-(m';j ean ) <mj/ 2 ' 3 <m„ 



3.(r(m^ em ), with m'™"" the mean mass of the three dark 
photon candidates and cr is the width of the signal as 
function of the dark mass which is taken from Monte 
Carlo (MC) simulation. The detection efficiency of 
Belle was modeled with MC simulations based on the 
GEANT4 package Il40il . The simulation includes all rel- 
evant properties of the sub-systems, including geomet- 
rical acceptance, charged particle identification, trig- 
ger efficiency, response of all detector modules, and 
selection criteria. The MC also includes information 
about inefficient individual detector modules. Belle can 
achieve, on average, a detection efficiency of 20 % and 
40 % respectively for 6 electron and 6 muon final states. 



5. Background estimation 

The background estimation is based on a data driven 
method. In this method, all combinations that have 
two pairs where the leptons are combined with their 
wrong-sign partner and one pair with opposite charge, 
(l~l~)(l + l + )(l + l~), are kept. The three masses are or- 
dered in decreasing order: ml > mf, > m^. The mass 
difference between the ml, the highest mass, and m\ 
the lowest mass: ml - m^ is then calculated. Fig- 
ure [T] shows the mass difference m\ - m\ as func- 
tion of the mass m l u for the 6 electrons (Figure [1} 
top) and 6 muons (Figure QJbottom) final states. For 
conciseness, the following notations are used: "same 
sign" for (l + 7~) and "opposite sign" for 
(l + r)(l + r)(l + r) (with / = e or fx). The scatter plot for 
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the opposite sign: (e + e~)(e + e~)(e + e~) - Figure QJtop- 
right and (jj. + fi~)(jj. + /j~)(fi + ijr) ~ Figure QJbottom-right 
have their signal region blinded (red bands), since a 
blind analysis technique is used and not all selection 
criteria have been validated. The signal region for the 
same sign scatter plots ((e + e~)(e + e + )(e~e~) - Figure Q} 
top-left and (/i + /U + )(yU~ - FigureQJbottom-left) 
is unblinded. 



data: 3e + 3e , m 1 - 1 .9 GeV/c 2 
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Figure 1 : versus m], — m\ for 6e final state (top) and 6/i final state (bot- 
tom). Top-left and bottom-left, opposite sign: (e + e~)(e + e~)(e + e~). Top-light 
and bottom-right, same sign: (e + e~)(e + e + )(e~ e~). The 2D histograms are di- 
vided into 20 slices, each slice is then projected on the mj, - m^-axis. The red 
airow (right) shows slice 10, corresponding to m), = 1.9 ± 0.1 GeV/c 2 and the 
direction of the projection. The red box on the left (top and bottom) corresponds 
to the blinded signal region. 



The background is then estimated for different ml mass 
regions as illustrated for the ml = 1.9 + 0.1 GeV/c 2 band 
by Figure [2] It is assumed that the backgrounds have 
the same shape in the side band region of the same sign 
distribution and the opposite sign distribution but not 
necessarily the same number of events. Therefore, the 
same sign distribution is normalized to the opposite sign 
distribution by a factor calculated for each m\ v The ex- 
pected background is then the scaled number of events 
counted in the signal region of the same sign distribu- 
tion. The signal obtained in the MC simulation is shown 
as a black curve in FigureQ](for the opposite sign). 



6. Expected sensitivity 

The expected background with no constraint on the 
impact parameters and the vertex positions is less than 
25 events and less than 5 events with a constraint on the 
impact parameters but no constraint on the vertex po- 
sitions for the six electron final states. No events are 
found in either case for the six muon final state. The 
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Figure 2: Projection on m l u - m 3 u for m\ t = 1.9 ± 0.1 GeV/c 2 . The "same 
sign' , -(e + e~)(e + e + )(e~e~) (black point) distribution has been normalized to the 
"opposite sign"- (e + e~)(e + e~)(eV~) (red point) distribution using the side band 
area. 
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Figure 3: Left, background estimation method verified successfully with MC 
Right, for experimental data, predicted BG is > 20 events for the electron final 
state. 



case with no constraint on the impact parameter and on 
the vertex position corresponds to displaced vertex po- 
sitions up to 80 cm from the interaction point and ac- 



cording to l2Ul24ll up to e ~ 10 6 . Here, e — — is the 
strength of the dark photon mixing with the Standard 
Model hypercharge gauge boson and a, the electromag- 
netic coupling. For the case with a constraint on the 
impact parameters, the decay length can go up to 6 mm 
i.e. up to e ~ 10~ 4 . 

Figure[3]-left shows a MC simulation test of the back- 
ground estimation method. The interaction e + e~ — » 
3e + 3e~ have been generated in phase space. The back- 
ground estimated from the same sign distribution is con- 
sistent with the number of events counted in the sig- 
nal region of the opposite sign for all niA candidates. 
Figure [fright shows the expected background deduced 
from the data same sign analysis. 

A statistical method based on the Feldman-Cousins 
approach 114 111 is used to then calculate preliminary up- 
per limits (90% CL) for a number of "observed" events 
equal to the estimated background events for the full 
luminosity and the prompt case. Figure [4] (left - six 
electron, right six muons final states) shows the prelim- 
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Figure 4: Sensitivity as function of the "Dark Photon" mass for different "Dark 
Higgs" mass hypothesis. Left: 6e. Right: 
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Figure 5: Sensitivity as function of the "Dark Higgs" mass for different "Dark 
Photon" mass hypothesis. Left: 6e. Right: 6fi. 



inary sensitivity limit for different "Dark Higgs" mass 
hypotheses. Figure [5] (left - six electron, right 6 muons 
final states) shows the preliminary sensitivity limit for 
different "Dark Photon" mass hypothesis. Due to the 
expected low level of background the sensitivity scales 
nearly linearly with the integrated luminosity. 



7. Conclusions 

The Dark Photon and the Dark Higgs are searched for 
in the mass ranges: 0.27 < niA < 3 GeV/c 2 and 
0.54 < m/,< < 10.86 GeV/c 2 . Based on control data sam- 
ples, it was found that the background is small, implying 
that the detection sensitivity scales nearly linearly with 
integrated luminosity. The expected preliminary Belle 
sensitivities have been shown. Results for data will be 
unblinded soon. 
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